from pylab import *

alpha = []

for i in range(1,101):
    a = loadtxt('alpha/alpha_'+str(i)+'.txt')
    tmp = []
    for j in range(20):
        tmp.append(a[j,0])
    alpha.append(tmp)
    
alpha = array(alpha)

for i in range(20):
    subplot(2,2,1+i/5)
    hist(alpha[:,i],bins=15,alpha=0.6,label=r'$\alpha_{'+str(i+1)+'}$ = '+str(round(mean(alpha[:,i]),4))+'$\pm$'+str(round(std(alpha[:,i]),4)) )
    xlim(-2.25,2.25)
    legend(loc='best',frameon=False)

show()
